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Abstract 



We analise the implications of the measurement of B and J/ip inclusive pt dis- 
tributions performed in pp collisions by the UA1 and CDF experiments. 



1 Introduction 



Heavy quark production in high energy hadronic collisions consitutes a fundamental arena 
for the study of perturbative QCD [|J]. The comparison of experimental data with the 
predictions of QCD provides a necessary check that the ingredients entering the evaluation 
of hadronic processes (partonic distribution functions and higher order corrections) are 
under control and can be used to evaluate the rates for more exotic phenomena or to 
extrapolate the calculations to even higher energies. The estimates of production rates 
for the elusive top quark rely on the understanding of heavy quark production properties 
within QCD. Following the initial encouraging agreement between the UAl measurements 
0, [| of inclusive bottom-quark and J /if) p t distributions and the best available QCD 
calculations ||, |^, ||, the subject has acquired a new particular interest after the latest 
measurements by UAl J7|, ^| and CDF || |TU|, |TTJ and after the recent studies of the small-x 



behaviour of hadro-production cross-sections [13, O, |14 



On the experimental side UAl has confirmed with the most recent analyses of the b 
sample the agreement between next-to-leading order (NLO) QCD and data, extending 
it to large values of pt|§. At the same time, however, UAl has found a significant 
disagreement between its latest J/ip Pt spectrum [|7| and the available leading order (LO) 
QCD calculations. Independently, CDF preliminary and published results indicate a clear 
discrepancy with respect to theoretical expectations for normalization and shape of the b 



and J ftp p t spectra || [10], [LI 



On the theoretical side, studies have indicated that large corrections to the NLO 



evaluation of the b cross section should be expected at very large energies [12, 13, 14]. A 



precise quantitative definition of how large these energies should be for these effects to 
become dominant is however still missing, and the question of whether these effects can 
indeed justify the findings by CDF is still open. 

In parallel, several measurements |T5L TE\ have improved the knowledge of the parton 
distribution functions (PDF) down to values of x of the order of 0.01, and the relative 
new sets of parametrizations have now become available W7[ ITS]. 

Both these latest developments could have some bearing on the measurements of the 
b and J/ip production rates at UAl and CDF, because the values of x probed by CDF 
are approximately a factor of 3 smaller than those explored by UAl. Therefore CDF is 
more sensitive than UAl to the possible uncertainties of the extrapolation to small- x of 
both PDF's and partonic cross- sect ions. 



Attempts have also been made [|19| |2(| to incorporate directly into the fits of gluon 
distributions the experimental information contained in the b cross section measurements 
by UAl and CDF. However when the new PDF measurements quoted above are included 
as an additional constraint, there seems not to be enough freedom to correct the theoretical 
prediction by the needed amount to produce a fully satisfactory agreement with CDF data. 

In this letter we will reconsider these data and use the most recent inputs to shed, if 
possible, some more light on these problems. 
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2 The inputs of the theoretical calculations 



We will start by discussing the ingredients of the calculation presented here and by com- 
paring them to the analyses quoted previously. 

J / ip production is evaluated by summing the result of inclusive b production followed 
by B — > J/ip + X decays, direct J/ip production and production via radiative decays 
of x states. These last two processes we will refer to as direct charmonium production, 
while the first one will be referred to as 6-decays. The matrix elements for charmonium 
production are taken from the original calculations of ref |H] , and therefore coincide with 
those used in previous studies from now on indicated by GMS). 

As an improvement with respect to GMS we will use: 

1. the full NLO matrix elements for the production of h quarks before their decay into 

2. the most recent inclusive B — > J/ip + X and B — > ip(2S) +X fragmentation spectra 

[0,113; 



3. the most recent parametrizations of PDF described in |T7j], exploring the dependence 
on the value of the 2- loop A4 within the one standard deviation range A4 = 215 ± 
60MeV. 

This third point, in particular, will cause a significant drop in the prediction of the absolute 
rates for charmonium compared to GMS, because the PDF's used in GMS had a much 
larger value of A4 (A4 _ ' oop = 400 MeV versus A4~' oop = 215 MeV as the central value of 
the new sets, corresponding to A4~' oop = 140 MeV.). 

As was done in GMS we will smear the 6-quark p t with a Peterson fragmentation 



function [23] before its decay into ?/>'s: 



dN z(l - z* 2 

Of 

dz [ez + (1 - z) 2 } 



« r^.\i 3212 . (2- 1 ) 



with e = 0.006. Furthermore we will assign to the 6-meson after fragmentation a mass of 
5.28 GeV, regardless of the input 6-quark mass. 

The central values we use for the 6-quark mass and the factorization/renormalization 
scale /i are m fe =4.75 GeV and /i 2 = /i 2 = Pt 2 + m b 2 Q We will consider in the following 
the effect of changing these parameters. The numerical values of other parameters used 
in the calculations are contained in Tables |l] and |2|. 

1 Notice that this prescription is slightly different from the one chosen in ||, where /1 2 = Ptmin + m b 2 = 
rriT 2 was used to estimate a(pt > pt mm )- This difference amounts to an effect of the order of —20% at 
the smallest values of pt . 
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2.1 Discussion of the uncertainties in the calculations 



Several factors contribute to the uncertainty of the calculations. We will discuss first 
those intrinsic to the perturbative approach and will consider later those related to ex- 
ternal inputs such as structure functions, Aq C d or potential models in the case of the 
charmonium production. 

First we consider b production. As is well known, the radiative corrections to b pro- 
duction are rather large and extremely unstable under changes in factorization scale ll, 
in particular at the larger CM energies. The standard way to establish how reliable the 
perturbative expansion is, is to vary the factorization scale ll within some range of the 
order of the hard energy scale relevant for the process considered. If the observed change 
in cross sections is large, one might think of selecting a best value for ll by fitting the 
measured rates. Since the dependence on ll varies as a function of the beam energy (indi- 
cating that the effect of yet higher order corrections is not the same at different energies) 
there is no reason a priori for which the expression for ll fitted at one value of the beam 
energy should be the same at different energies. 

Therefore the uncertainty related to the choice of ll cannot be removed by performing 
independent measurements at different energies. 

Likewise, there is no guidance on what is the proper range within which to allow ll 
to vary. If [lq is the typical energy scale of a given process (say the transverse mass of a 
heavy quark or the E t of a jet), it is customary to vary llq/2 < 11 <2lio- There are several 
indications, however, that when working at a fixed order in perturbation theory scales 
significantly smaller than the natural scale are needed to reproduce the data. As examples, 
we quote the cone-size dependence of the jet E t distributions in hadronic collisions [pj 



or the jet multiplicity distributions in e + e~ [EJj]. Wherever all-order calculations have 
become available, these indicate that the resummation of leading and sub-leading large 
logarithmic terms at any order in the perturbative expansion restores the insensitivity to 
li and allows ll to be chosen of the order of the natural scale llq . 



We therefore believe it is legitimate to push ll to values as low as possible, compatibly 
with the range allowed by the PDF parametrizations. In our case, we will consider the 
range m-r/4 < ll <itlt, which will give us ll 2 > 5 GeV 2 for all values of p t probed by the 
CDF data, and therefore does not include regions of Q 2 which are not under the control 
of the DIS data. 

Similar considerations apply to the case of direct quarkonium production of J/ip's. 
Here the situation is even worse, because only the LO production processes are available 
and the expected ll dependence is even more significant. In GMS the effect of possible 
higher order terms was parametrized in terms of a constant K factor, chosen to have a 
value of 2 to reproduce ISR data. As mentioned above, however, there is no reason a 
priori why this same value for K should apply at the significantly higher energies used 
at UA1 and CDF. Once again, therefore, we will choose to probe the possible effects of 
higher order terms by selecting different values of ll. In accordance with the choice made 
for the b production, we will choose < ll <m>T, where in this case my represents 

the transverse mass of the quarkonium state. As will be shown later on, different values 
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of fi will not only change the absolute normalization but will also affect the shape of the 
p t distribution of J/^'s. 

In addition to the above uncertainties, one should add the uncertainty in the evaluation 
of the parameters of the quarkonium states entering the estimate of their production cross 
section. As an example, we quote a recent study |27| of |i?'(0)| 2 - the first derivative of 
the wave function at the origin for P-wave states - indicating a value for Xc states which 
is approximately 50% higher than what obtained from previous models p8 |. The values 
used in this paper are given in Table 1, and follow the old results of pH . 



2.2 Numerical Results 

We collect the results obtained for the b and J/ip cross sections as a function of the 
various input parameters in a series of tables. As a standard reference we will use sets DO 
and D of the recent MRS PDF parametrization | I7H . We will use two values of Aqcd-, 



A°=215 MeV and A^=275 MeV, corresponding respectively to the central value and to 
one standard deviation above the central value obtained from the fit. Tables [3] and [| 
contain the bottom quark p t distribution integrated above a given p t at 1.8 TeV and for 
the two extreme values of /z, /i=mT and /i=mr/4. The quark is required to satisfy \y\ < 1, 
to allow comparison with the CDF data. Tables [| and [7] contain the same information, 
but at 630 GeV and with \y\ < 1.5, to allow comparison with UA1 data. 

Several comments are in order. First of all notice that while the use of the more 
singular set of structure functions leads to larger values of the total cross sections (p t > 0) 
at 1.8 TeV, the opposite happens at 630 GeV. This is because the higher density of gluons 
at small x described by set D forces via momentum sum rules a depletion at larger values 
of x. Since UAl is sensitive to larger values of x, an overall decrease in the total cross 
section is observed. 

Notice also, on the other hand, that even at CDF the singular gluon parametrization 
D- will give a cross section smaller than the set DO as soon as we consider transverse 
momenta of the b above 10 GeV - which is the region where most of the CDF data are. 
Since above 10 GeV the shapes of the integrated pt distributions for the two parametriza- 
tion DO and D are similar, this indicates that the measurement of the cross section for b 
production in this region cannot be reliably used to extract via extrapolation to pt=0 the 
total b production cross section. For example, while the region p t > 10 GeV represents 
10% of the total cross section according to DO and using Aqc,d=215 MeV, the same re- 
gion represents only 7% of the total according to D and using Aqcd=275 MeV. A similar 
exercise at the UAl energy indicates a more reliable extrapolation. 

As already indicated in ||, the dependence on the value of the b mass is not significant. 
In Table |5| we show a comparison between the integrated b pt distribution obtained using 
m&=4.5 and m;,=4.75 GeV. The difference is of the order of 20% for the total cross section, 
but becomes negligible for p t > 10 GeV. 

In Table § we present the integrated p t distribution of J/ip mesons, calculated at 
CDF energy and divided into the direct quarkonium and B decay contributions. The 
relative fraction due to B decays, indicated by Jb, is also shown as a function of the pt 
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threshold, and the dependence on Aqcd and \i is studied by considering the central case of 
Aq C - £ )=215 MeV, n=m T and the extreme case of Aqc<£>=275 MeV, n=m T /A. A priori there 
is no reason why the same factorization scale should be used for the two contributions, as 
the two physics processes are entirely different. Furthermore the B decay is evaluated at 
NLO, while as mentioned previously only the LO terms are available and included in the 
quarkonium term. Nevertheless we take here the value of /i for the two processes to be 
the same, in order to extract and indicative range of values for The value of plays 
an important role in the experimental determination of the B cross section out of the 
measurement of the inclusive J/ip rate, and the range of values exhibited by the tables 
indicates what is the systematic uncertainty that one should expect in deriving fs from 
the theory. Parallel results for UAl energy are shown in Table || 

The most important thing to notice about these tables is the fact that the B contri- 
bution only changes within a factor of 2 by changing /i, while a variation ranging from 
a factor of 7 to 10, depending on p t , is observed for the charmonium case. This indi- 
cates that the LO prediction for direct charmonium is very poor, and very large NLO 
corrections should be expected. 

3 Comparison with the data and discussion 

We will now compare UAl and CDF data with the results obtained so far. For these 
comparisons we use the results obtained with the DO PDF set and with the two different 
choices (/i, Aqcd) — (wr, 21BMeV) and (h,Aqcd) — (m T /4:,275MeV), which provide an 
acceptable upper and lower limit to the band of current theoretical uncertainty relative 
to the inputs discussed above. 

Figs. [l| and |2| show the UAl measurement of the inclusive 6-quark and B meson p t 
distribution integrated above a given threshold p™ in and with \y\,,B\ < 1-5. The agreement 
between data and theory observed in ref.|| is confirmed, even though the central value 
of the prediction has dropped by almost 50% as a consequence of the smaller value of A 4 
in the central MRS fit compared to the central DFLM fit (215 versus 260 MeV). Also the 
S-meson spectrum is well consistent with what expected from a Peterson fragmentation 
model, as anticipated in ||. Notice, however, that a priori there is no guarantee that the 
Peterson model should work for values of p t of the order of the B mass, as in this region 
corrections to factorization could be significant. 

Fig.^ shows the inclusive p t differential distribution for J/i/j's produced with p t > 5 
GeV and \y\ < 2. We superimpose the contributions from direct charmonium production, 
b-decays and the sum of the two. The data fall all inside the theoretical band. Since as 
mentioned above there is no reason to expect the expression for /i to be the same for the 
two contributions - while for the sake of simplicity we imposed this in adding the separate 
terms in the figure - a better fit to the data could be obtained by choosing /i=m^/4 for 
the charmonium and n=rriT for the B production. 

Fig. [| shows the integrated p t distribution of b quarks with \y\ < 1 from CDF, com- 
pared to the results of the NLO calculation. The data are taken from published results 
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as well as from recent public presentations ||, |10|, [TT|. As already observed in || there 
is a clear excess in the observed rate at small p t . At larger values of p t , in the region of 
the inclusive b — > / + X measurements, the data are consistent with the upper estreme of 
the theoretical band. A similar feature is observed in the tp differential pt measurement, 
shown in fig. [| 

Equally worrysome is the comparison between theory and data in the case of the pt 
spectrum of the ip(2S), shown in fig. H As noted in |10|, the expected contribution 
from direct quarkonium production is heavily suppressed. We confirm this estimate, and 
verified that it remains true even allowing for the variation of /x within the /io/4< \i </U 
range. 

Is it possible to explain the patterns observed by UAl and CDF in a unified fashion 
by invoking generic small-x effects, either from PDF's or from violation of factorization? 
Rather than studying this question by directly attempting to modify the gluon densities, 
as done in refs. [TP], |2(J , we will address it here by considering the following quantity: 

a{x g < x;p t b > p t mm ) = [ dx g da(Ptb > Vrn) , (3.1) 

namely the contribution to the integrated pt distribution coming from partons with mo- 
mentum fraction smaller than a given value of x. We plot this variable as a function 
of x and for different values of p t mm (b) in Fig. [7] for UAl and CDF. We only integrated 
over b quarks within the regions of acceptance of the experiments, namely \yb\ < 1.5 for 
UAl and \yb\ < 1 for CDF. Since the contribution to the cross-sections due to the qq and 
qg initial states are negligible for the relevant regions of p t we are concerned with, we 
limited ourselves to the gg process and normalized the curves to the value of 1 at x — 1. 
Therefore the plotted functions represent the fraction of cross-section due to gluons with 

Xg < X. 

The first thing to notice is that the distribution corresponding to p t > 5 GeV at UAl 
lies between the curves for p t > 10 and p t > 20 GeV at CDF, consistently with the 
factor of 3 difference in beam energy. The second thing to notice is that at CDF energies 
the contribution to the cross section for p t > 10 GeV from the region x < 0.01 is less 
than 20%. Furthermore no contribution at all comes from the region x < 0.003. We 
verified that different fits of the NMC and CCFR data, obtained in Ref. [|18||, give gluon 



densities which differ, over the relevant kinematic range, by no more than 10% from the 
MRSD0 set used here. Since all of these gluon parametrizations do not differ significantly 
from previous extrapolations, we conclude that the knowledge of the gluon density in the 
relevant region 0.1 > x > 0.01 and Q > 5 is today rather solid. We therefore expect that 
only dramatic changes in the guon densities in the region 0.003 < x < 0.01 will lead to a 
change of a factor of 2 in the cross section integrated above Pi=10 GeV. 

Therefore while it is tempting to conjecture that the ignorance about the behaviour 
of the gluon densities at small-x could explain the discrepancy between the overall rates 
measured by UAl and CDF and the difference in slope of the CDF spectra compared to 
theory, we find no evidence that this assumption is justified. Rather, we find that the 
region x g < 0.01 is marginal in the production of b quarks or ip's passing the required 
acceptance and p t cuts imposed by the two experiments. 
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The effect of the small-x corrections to the partonic cross-section considered in |14 
is more difficult to estimate. In fact these phenomena alter the kinematic connection 
between p t and x, since they predict that initial state gluons with a given momentum 
fraction x can have a pt non negligible w.r.t. xEb eam . This is equivalent to having an 
intrinsic pt of the order of the scale of the hard process itself, namely rrib- As a result, the 
region with x g < 0.01 could provide a significant contribution to the rate for p\ > 10 GeV, 
thanks to the transverse momentum smearing induced by this sort of small-x primordial 
p t . Even though it was found in ref. [0 that these small-x effects can add at most 50% 



of the NLO contribution to the total b cross section at 1.8 TeV, no explicit indication 
is given on the ^distribution of this additional 50%. Since the cross section observed 
experimentally (p tb > 8.5 GeV) represents of the order of 10% of the total rate at NLO, 
we cannot exclude that the p t smearing induced by these effects be responsible for the 
factor of 2-3 discrepancy observed between data and NLO predictions. Notice that the 
hypotesis of a p t smearing would help understanding not just the rate deficiency, but 
also the apparent difference in shape between NLO and data. A quantitative statement 
regarding these possiblities will only come from more explicit studies along the lines of 
ref.P|. 



While we await for more explicit calculations, it might be worth exploring some ad- 
ditional consequences of this scenario. In addition to trying to push the measurement of 
the b cross section to even smaller values of p t , it would be important to study correla- 
tions between the pair of b quarks. NLO calculations exist for these correlations [pSfl . If 



the small-x effects were to behave as indicated previously, we would expect to observe a 
flattening of the A0 and pt h distributions w.r.t the NLO prediction. Here A<f) represents 
the difference in azimuth between the b and the b, and p t bb represents the transverse mo- 
mentum of the pair. The flattening would be caused by the additional intrinsic p t due to 
the gluon transverse momentum. 



Measurements of the A0 correlations have been performed by UAl [[3(J, indicating a 



good agreement with the NLO calculation [29]. This result does not resolve the issue, 



though, because the agreement of the NLO b cross section with the data suggests that 
the energy at UAl is below the threshold for the possible onset of these new small-x 
phenomena. 



4 Conclusions 

After allowing for rather generous estimates of the theoretical uncertainties involved in 
the calculations currently available for b and J/ip production in hadronic collisions, we 
conclude that the most worrisome points of discrepancy can be summarised as follows: 

1. The production of direct charmonium both at CDF and UAl is much more abundant 
than would be obtained from the LO calculation using a standard value of /i=m*r. 
The J/ip can be explained by using pbit/4, which however gives a rate 8-10 times 
larger than for fi=rriT, indicating a rather unstable perturbative expansion. However 
this is not sufficient to explain the rate of ip{2S) production. 
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2. b production at CDF for values of pt around 10 GeV is significantly larger than 
can be accomodated by current estimates of the higher order effects or by possible 
structure function uncertainties. Using the extreme value of /x=m T /4 is not sufficient 
to explain all the data points, the discrepancy being still larger than a factor of 2. 
We cannot however exclude that the solution be in the large smearing induced by a 
small-x intrinsic pt of the initial state gluons. At larger values of pt we believe that 
the consistency between data and theory is acceptable. 



What other effects could be responsible for the remaining discrepancies? It should be 
noticed that the two points above might not be uncorrelated. In fact the absolute nor- 
malization of the two CDF points at lower p t , coming from the measurement of inclusive 



J/ip and ip(2S) rates [1C], relies on two assumptions: (i) that all of the ip(2S) come from 
B decays, and (ii) that the J / ip fraction f B is known. The right hand side of Table § - 
which represents the choice of parameters which comes closer to representing the CDF 
J/ij} spectrum - suggests a value for fs which is significantly smaller than the central 
value used by CDF (namely 37% vs. 63± 17% for p t (ip) > 6 GeV Q). This would 
decrease the effective b cross section by a factor of 50%. In addition, the new processes 
responsible for the large K factor apparent in J/ijj production might affect ip(2S) produc- 
tion and could provide enough rate to reduce the b rate extracted from the assumption 
that B decays are the only source of ip(2S). 

Notice that /# could in principle be extracted experimentally, for example by sepa- 
rating the direct J/ip's from those due to B decays via the observation of the displaced 
vertex from which the ip orginates - due to the long B lifetime. UAl measured fs by 
assuming that direct J/ip's are isolated while J/^'s from B decays are not, and studying 
the isolation of the J/^'s in the data. This assumption however might not be correct if 
other production mechanisms were responsible for direct quarkonium production, such as 
for example gluon — > J/tp fragmentation [^TJ . 

It is very reasonable to expect that at some value of pt the dominant production mech- 
anism for charmonium states will indeed be via gluon fragmentation. The main reason 
being that direct production as described by the LO mechanisms inhibits production at 
large p t via a form factor suppression (the probability that a charmonium bound state 
will hold together when produced directly in an interaction with a large virtuality scale is 
highly suppressed). The fragmentation functions for the creation of S"-wave charmonium 



(r) c and J / ip) in a gluon shower have recently been calculated [|3T| and work on the creation 



of P-wave states (x) is in progress |32|]. It will be interesting to use these calculations in 



order to extract the fragmentation contribution to charmonium production in the regions 
of pt explored experimentally, and verify whether these new processes can account at least 
in part for the large observed K factor. The experimental detection of non-isolated J/^'s 
from a primary vertex - and therefore presumably not coming from B decays - could 
provide a strong indication that these processes are indeed present. 

Similar measurements of the decay-vertex position of the ip{2S) would provide evidence 
in favour or against the current belief that most of them come from B decays. Once again 
the gluon fragmentation contribution to production of this charmonium state could turn 
out to be significant, and would manifest itself with a signal of non-isolated prompt ip(2S). 
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Similarly interesting would be a separate measurement of the x Pt spectrum, which 
is expected to be dominated by direct production rather than B decays. A preliminary 
measurement by CDF |TTJ reports BR{jp — > /i + /i~)x a(xc ~^ i?Ti Vt x > 7GeV; \rj\ < 0.5) 
= 3.2 ± 0.3 ± 1.2 nb. Both xi an d X2 are here included. This can be compared with 
the range O.QAnb < a < 5.1nb obtained using the calculation described in the previ- 
ous sections and the two extreme choices (/1, Aqct>) = (%, 215 MeV) and (fi, Aqcd) = 
(mr/4, 275MeV). Using the above cross section and using the inclusive B — > Xci branch- 
ing ratio of 0.54 ± 0.21% |3"3" |, we estimate that only a fraction of the order of 10% or less 
- depending on p t - of the x' s come from B decays. Since the production mechanisms for 



Xi and X2 are different even at LO [^TJ, a separate measurement of the two states would 
be welcome, even though their closeness in mass makes it very hard to separate one from 
the other in practice. 

It would also be interesting to evaluate the effects of resumming some of the leading 
and next-to-leading corrections to the evolution of the initial state, using the calculation 
of the gg — > color-singlet NLO form factor calculated in ref. 133 



An experimental measurement of the production cross section and p t spectrum for 
T states would be very useful in understanding the quarkonium production mechanisms 



35fl . In this case, in fact, one would have at least three advantages: (i) the masses 



involved are larger and presumably both the non-relativistic approximation involved in 
the determination of the quarkonium wave function and the QCD perturbative expansion 
would work much more reliably than for charmonium; (ii) the signal does not have a 
contamination similar to the one due to B decays; (iii) the p t spectrum could hopefully 
be extended to very small values of pt, possibly even to pt=0, thanks to the large mass of 
the T and the rather large momentum and easier detection of the decay muons. 
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fe(0) | 2 


|i?' 1P (0)| 2 


0.7 


0.4 


0.006 



Table 1: Values of wave functions used in the evaluation of the charmonium cross 
sections. 



BR(B -> ipX) x BR(i[> -> 
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Pt 


MRSDO 


MRSD- 


[\je v ) 


A 

A 


A 


A 

A 


A 


n 
U 


1 1 j|p i Oil 


1 OAT?_i_CiA 
1.04ri/+U4 


1.00I1/+U4 


1.0 ( I1/+U4 


r 



4.0UJI/+U0 


O.ZZil/+Uo 


4.0UJI/+U0 


O.z / J1/+U0 


1 n 
1U 


i.uon/-ruo 






± .uyj_/-rUo 


10 


3.16E+02 


3.56E+02 


2.70E+02 


3.11E+02 


20 


1.15E+02 


1.32E+02 


9.69E+01 


1.12E+02 


25 


4.97E+01 


5.58E+01 


4.20E+01 


4.75E+01 


30 


2.40E+01 


2.72E+01 


2.06E+01 


2.34E+01 


40 


6.94E+00 


7.85E+00 


6.05E+00 


6.80E+00 


50 


2.60E+00 


2.79E+00 


2.16E+00 


2.38E+00 


59 


1.21E+00 


1.30E+00 


1.03E+00 


1.14E+00 



Table 3: Integrated bottom quark p t distribution at 1.8 TeV. m b =4.75 GeV, /i=/x , 
A = 215 MeV, A+ = 275 MeV. 



Pt 


MRSDO 


MRSD- 


(GeV) 


A 


A+ 


Ao 


A + 





2.17E+04 


3.03E+04 


2.69E+04 


3.71E+04 


5 


7.23E+03 


9.39E+03 


7.68E+03 


1.00E+04 


10 


1.83E+03 


2.27E+03 


1.66E+03 


2.09E+03 


15 


5.78E+02 


7.11E+02 


4.93E+02 


6.10E+02 


20 


2.23E+02 


2.68E+02 


1.83E+02 


2.25E+02 


25 


9.86E+01 


1.17E+02 


7.98E+01 


9.75E+01 


30 


4.84E+01 


5.63E+01 


3.89E+01 


4.72E+01 


40 


1.44E+01 


1.65E+01 


1.23E+01 


1.38E+01 


50 


5.28E+00 


6.04E+00 


4.48E+00 


4.87E+00 


59 


2.23E+00 


2.49E+00 


2.04E+00 


2.27E+00 



Table 4: Integrated bottom quark p t distribution at 1.8 TeV. mf,=4.75 GeV, /z=/zo/4, 
A = 215 MeV, A + = 275 MeV. 
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Pt 


Ao 


A+ 


(GeV) 


m b =4.5 GeV 


m fc =4.75 GeV 


m 5 =4.5 GeV 


m 6 =4.75 GeV 





2.6E+04 


2.2E+04 


3.8E+04 


3.0E+04 


5 


8.0E+03 


7.2E+03 


1.0E+04 


9.4E+03 


10 


1.9E+03 


1.8E+03 


2.4E+03 


2.3E+03 


20 


2.3E+02 


2.2E+02 


2.8E+02 


2.7E+02 



Table 5: Mass dependence of the integrated bottom quark p t distribution at 1.8 
TeV. MRSDO parton distributions, ^=A*o/4, A = 215 MeV, A + = 275 MeV. 



Pt 


MRSDO 


MRSD- 


(GeV) 


A 


A + 


A 


A + 





6.016E+03 


7.110E+03 


5.584E+03 


6.606E+03 


5 


1.838E+03 


2.149E+03 


1.582E+03 


1.848E+03 


10 


2.953E+02 


3.390E+02 


2.482E+02 


2.857E+02 


15 


6.353E+01 


7.183E+01 


5.378E+01 


6.188E+01 


20 


1.776E+01 


2.009E+01 


1.566E+01 


1.778E+01 


25 


6.148E+00 


6.714E+00 


5.271E+00 


6.110E+00 


30 


2.419E+00 


2.620E+00 


2.143E+00 


2.378E+00 


40 


4.654E-01 


5.086E-01 


4.616E-01 


5.098E-01 


50 


1.312E-01 


1.480E-01 


1.324E-01 


1.381E-01 


59 


4.672E-02 


5.346E-02 


4.242E-02 


4.552E-02 



Table 6: Integrated bottom quark p t distribution at UA1. m b =4.75 GeV, fJ,=Ho, 
A = 215 MeV, A + = 275 MeV. 
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Pt 


MRSDO 


MRSD- 


(VoreV ) 


A 


A 

A + 


A 


A 


U 






1 on ir i ri/i 


1 CX A *~7T? 1 n A 

1.04 ^Hj+U4 


rr 



o a i niT 1 i no 

3.619E+03 


/i h 7nci? 1 no 

4.705E+03 


nr'nTT 1 1 no 

3.069E+03 


no txTT 1 1 no 

3.985E+03 


1 n 

10 


/■ / 1 /" T7 1 i r^o 

6.036E+02 


t /i Tnir 1 1 no 

7.479E+02 


A no/^T? 1 no 

4.986E+02 


f* oaot? 1 no 

6.202E+02 


15 


i qi nx? i no 


1.0UU1I/+UZ 


1 1 qqtt i no 




on 


3.635E+01 


4.344E+01 


3.193E+01 


3.666E+01 


25 


1.184E+01 


1.391E+01 


1.092E+01 


1.235E+01 


30 


4.452E+00 


5.139E+00 


4.210E+00 


4.575E+00 


40 


8.034E-01 


9.621E-01 


7.467E-01 


8.312E-01 


50 


1.723E-01 


1.828E-01 


1.765E-01 


1.579E-01 


59 


4.111E-02 


4.241E-02 


4.878E-02 


4.510E-02 



Table 7: Integrated bottom quark p t distribution at 630 GeV. m b =4.75 GeV, 
/i=/i /4, A = 215 MeV, A+ = 275 MeV. 



min,ib 
Pt 


A , /i=/i 




A 


+, n=Ho/± 




(GeV) 


a B - BR (nb) 


a x - BR (nb) 


fs (%) 


Ob- BR (nb) 


o x - BR (nb) 


fB (%) 


3 


2.6E+00 


5.4E+00 


32 


5.4E+00 


3.9E+01 


12 


4 


1.7E+00 


1.9E+00 


46 


3.4E+00 


1.5E+01 


19 


5 


1.1E+00 


7.6E-01 


58 


2.2E+00 


6.0E+00 


26 


6 


6.7E-01 


3.4E-01 


66 


1.4E+00 


2.8E+00 


33 


8 


2.9E-01 


8.7E-02 


77 


6.3E-01 


7.7E-01 


45 


10 


1.4E-01 


2.9E-02 


83 


3.2E-01 


2.6E-01 


55 


12 


7.8E-02 


1.2E-02 


86 


1.7E-01 


1.0E-01 


62 


14 


4.4E-02 


4.9E-03 


89 


1.0E-01 


5.0E-02 


67 


16 


2.6E-02 


2.5E-03 


91 


6.1E-02 


2.5E-02 


71 


18 


1.7E-02 


1.3E-03 


92 


3.9E-02 


1.3E-02 


75 


20 


1.1E-02 


7.3E-04 


93 


2.5E-02 


6.9E-03 


78 


25 


4.2E-03 


1.9E-04 


95 


9.4E-03 


1.9E-03 


82 


30 


1.9E-03 


4.9E-05 


97 


4.1E-03 


5.2E-04 


88 



Table 8: Integrated ip p t distribution from B decays, from charmonium production 
(X + V>) and relative B fraction at 1.8 TeV. MRSDO, A = 215 MeV, A + = 275 MeV. 
BR(J/^> — > /i + /i~) included. 
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min.tb 
Pt 


Aq , /!— /iQ 




A 


+, /x=/i /4 




(GeV) 


o B - BR (nb) 


o x - BR (nb) 


fs (%) 


o B - BR (nb) 


o x - BR (nb) 


fs (%) 


3 


3.0E+00 


3.4E+00 


46 


7.8E+00 


2.9E+01 


21 


4 


1.6E+00 


1.2E+00 


56 


4.0E+00 


1.1E+01 


27 


5 


8.3E-01 


4.7E-01 


63 


2.1E+00 


4.3E+00 


33 


6 


4.6E-01 


2.1E-01 


69 


1.2E+00 


1.9E+00 


37 


8 


1.6E-01 


5.1E-02 


75 


4.0E-01 


5.1E-01 


43 


10 


6.1E-02 


1.6E-02 


79 


1.6E-01 


1.7E-01 


47 


12 


2.7E-02 


6.2E-03 


81 


6.7E-02 


6.7E-02 


50 


14 


1.3E-02 


2.7E-03 


83 


3.1E-02 


2.9E-02 


51 


16 


6.6E-03 


1.2E-03 


84 


1.6E-02 


1.4E-02 


53 


18 


3.6E-03 


6.1E-04 


85 


8.4E-03 


7.2E-03 


53 


20 


2.0E-03 


3.2E-04 


86 


4.5E-03 


3.9E-03 


53 


25 


5.8E-04 


7.2E-05 


88 


1.2E-03 


9.2E-04 


55 


30 


2.0E-04 


1.8E-05 


92 


3.6E-04 


2.3E-04 


61 



Table 9: Integrated ip p t distribution from B decays, from charmonium production 
(X + VO and relative B fraction at 630 GeV. MRSD0, A = 215 MeV, A+ = 275 MeV. 
BR(J/ip — > /i + /U~) included. 
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Figure 1: Integrated b p t distribution at 630 GeV: UA1 data f§ versus NLO QCD. The 
J/ip point assumes a B fraction in the inclusive J/ip sample of 31±12% J/J. 



Figure 2: Integrated B meson p t distribution at 630 GeV: UA1 data versus NLO QCD. 



Figure 3: Differential J/ip pt distribution at 630 GeV: UAl data versus different QCD 
contributions, as shown in the legend. 



Figure 4: Integrated b p t distribution at 1.8 TeV: CDF data [11 versus NLO QCD. The 



J/ip point assumes a B fraction in the inclusive J/ip sample of 63±17% |i 1 



Figure 5: Differential J/ip p t distribution at 1.8 TeV: CDF data versus different QCD 
contributions, as shown in the legend. 



Figure 6: Differential ip(2S) p t distribution at 1.8 TeV: CDF data versus total QCD. The 
different contributions from direct production and B decays are labeled as in the previous 
Figure. 



Figure 7: Contribution to the b cross section above given p t thresholds as a function of 
the gluon momentum fraction x. 
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